Resonantly driven wobbling kinks 
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The amplitude of oscillations of the freely wobbling kink in the (f> 4 theory decays due to the 
emission of second-harmonic radiation. We study the compensation of these radiation losses (as 
well as additional dissipative losses) by the resonant driving of the kink. We consider both direct 
and parametric driving at a range of resonance frequencies. In each case, we derive the amplitude 
equations which describe the evolution of the amplitude of the wobbling and the kink's velocity. 
These equations predict multistability and hysteretic transitions in the wobbling amplitude for each 
driving frequency — the conclusion verified by numerical simulations of the full partial differential 
equation. We show that the strongest parametric resonance occurs when the driving frequency 
equals the natural wobbling frequency and not double that value. For direct driving, the strongest 
resonance is at half the natural frequency, but there is also a weaker resonance when the driving 
frequency equals the natural wobbling frequency itself. We show that this resonance is accompanied 
by translational motion of the kink. 

PACS numbers: 05.45.Yv 
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I. INTRODUCTION 



The kink of the <f) A equation has a mode of internal 
oscillation, commonly referred to as the wobbling mode. 
To check whether the resonant excitation of the wob- 
bling mode could provide a channel for pumping energy 
into a kink-bearing system, several authors have stud- 
ied the dynamics of (f) A kinks subjected to resonant di- 
rect or parametric driving and damping [H, 0, HI 5 • The 
damped-driven (f> A theory arises in a variety of physical 
contexts, in particular in the description of topological- 
defect dynamics in media with temporally [f| and spa- 
tially [6[ modulated parameters or in the presence of 
fluctuations Q. Examples include the drift of domain 
walls in ferromagnets in oscillatory magnetic fields Q; 
the Brownian motion of string-like objects on a periodi- 
cally modulated bistable substrate [lj3] ; ratchet dynamics 
of kinks in a lattice of point-like inhomogeneities [ll| and 
rectification in Josephson junctions and optical lattices 
.12). The damped (j> A equation driven by noise was used 
to study the production of topological defects during the 
symmetry-breaking phase transition (l3j and spatiotem- 
poral stochastic resonance in a chain of bistable elements 

The mathematical analysis of the damped-driven kinks 
started with the work of Kivshar, Sanchez and Vazquez 
Q who studied the discrete parametrically pumped (j> 4 
system. Assuming that the driving frequency lies above 
the phonon band and using the method of averaging, they 
have discovered the phenomenon of kink death for suf- 
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ficiently large driving strengths. Next, in an influential 
paper |8y , Sukstanskii and Primak considered the contin- 
uous </> 4 equation with a combination of direct and para- 
metric driving. (See also a related discussion in Q.) Us- 
ing a variant of the Lindstedt-Poincare technique where 
the velocity of the kink is adjusted so as to eliminate sec- 
ular terms at the lowest orders of the perturbation expan- 
sion in powers of the wobbling amplitude, they detected a 
slow unidirectional motion of the kink. In their analysis, 
Sukstanskii and Primak were not concerned with terms 
higher than quadratic in the amplitude of the wobbling; 
in fact, their approach is not suitable to deal with secular 
terms at the e 3 -order of the expansion. Neither can it be 
utilised in the case of the resonant driving frequency. 

The resonant situation was studied by Quintero, 
Sanchez and M ertens 0, 0, 0] who employed the method 
of projections. The method assumes a specific functional 
dependence of the kink on the "collective co-ordinates" , 
which in this case are the width and position of the 
kink. Inserting the chosen ansatz in the partial differ- 
ential equation and projecting the result onto the neu- 
tral modes associated with the two degrees of freedom, 
one obtains a two-dimensional dynamical system, a sim- 
plification from the infinite degrees of freedom present 
in the original partial differential equation. (In the un- 
damped undriven situation, the collective coordinate ap- 
proach was pioneered by Rice and Mele [H, Il6j|.) The 
major advantage of the method is that if the collective 
coordinates have been chosen such that they capture the 
essentials of the dynamics, one should be able to uncover 
the very mechanism of the observed nonlinear phenom- 
ena — which is otherwise concealed by an infinite number 
of degrees of freedom (see e.g. 17:]). The drawback of 
the collective coordinate approach, however, is that one 
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essential and which can be omitted without a qualitative 
impact on the kink's dynamics. Specifically, the role of 
radiation (which is neglected in the approach in ques- 
tion) is not obvious and therefore one has no a priori 
guarantee that a radiationless ansatz is adequate. An- 
other disadvantage is that the resulting two-dimensional 
dynamical system is amenable to analytical study only 
in a very special case (when the damping is set to zero 

In this paper, we approach the resonantly driven wob- 
bling kink from a complementary perspective. Instead 
of trying to guess the most pertinent set of collective 
variables, we construct the wobbler as a singular per- 
turbation expansion using a sequence of space and time 
scales. The nonsecularity conditions yield equations for 
the slow-time evolution of the wobbling amplitude and 
kink velocity. This asymptotic procedure has already 
been implemented for the unperturbed c/> 4 equation [21] : 
here we extend it to include damping and driving terms. 
We consider both external (direct) and parametric driv- 
ing, at several resonant frequencies. This includes the 
cases considered by Quintero, Sanchez and Mertens in 
@) H, 13] • Although the multiscale expansion cannot crys- 
tallise the "minimum set" of degrees of freedom account- 
able for the observed behaviour, it does not suffer from 
the arbitrariness associated with the choice of collective- 
coordinate ansatz. The multiscale expansion is asymp- 
totic, i.e. only valid for small amplitudes of the wobbling 
mode; however, in the small-amplitude limit the expan- 
sion provides a faithful description of the wobbler, in- 
dependent of any assumptions and mode pre-selections. 
Importantly, it does not neglect the radiation. 

In all four cases of the direct and parametric driv- 
ing considered in this paper, we derive an autonomous 
system of equations for the amplitude of the wobbling 
mode and velocity of the kink. In each of the four cases 
this dynamical system turns out to exhibit stable fixed 
points corresponding to the nondecaying wobbling of the 
kink. In some parameter regimes, the amplitude of sta- 
ble wobbling is nonunique and may undergo hysteretic 
transitions between two nonzero values. In another case, 
the wobbling is necessarily accompanied by translational 
motion of the kink. The conclusions of our asymptotic 
analysis have been verified in direct numerical simula- 
tions of the corresponding partial differential equation. 
The numerical procedure we shall be using throughout 
this paper was specified in [2l] , 

An outline of this paper is as follows. In sections HTl and 
IIIII we study the parametrically driven wobbling kink. In 
section [TT] the frequency of the driver is chosen near the 
natural wobbling frequency of the kink while in IIIII we 
take the forcing frequency close to double that value. In 
the next section (IIV|) . we compare the mechanisms that 
are at work in each of the two cases. Subsequently, we 
consider the kink driven directly — first near half of its 
natural wobbling frequency (section [V| and then close to 
the wobbling frequency itself (section IVI[) . Our conclu- 
sions are summarised in section lYlII Here, in particular, 



we rank the four resonances according to the amplitude 
of the resulting stationary wobbling and according to the 
width of the resonant frequency range. 

II. 1:1 PARAMETRIC RESONANCE 

A. Asymptotic Multiscale Expansion 

We start with the parametric driving of the form 
\4>tt - +~/&-[l + h cos(nt)] </> + 3 = 0. (1) 

This type of a driver was previously considered by Quin- 
tero, Sanchez and Mertens [ij. The driving frequency 
f2 is assumed to be slightly detuned from ojq, the linear 
wobbling frequency of the undriven kink: 

n = qjq(i + p). 

We remind the reader that ujq — \/3 [21]. Introducing a 
small parameter e (which will be used to measure the am- 
plitude of the wobbling mode in what follows) we choose 
the detuning in the form 

p = t 2 R, 

where R is of order one. Since the frequency of the free 
nonlinear wobbling is smaller than wo {see Eq.(44) in 
[2l]}. we expect that the strongest resonance will occur 
not when f2 = ojq but for some small negative p. [This 
will indeed be the case; see Eq. (l29l) below.] 

Next, the parameters 7 > and h > are the small 
damping coefficient and driving strength, respectively. 
We choose the following scaling laws for these param- 
eters: 

7 - e 2 r, h = e 3 H, (2) 

where T and H are quantities of order 1. This choice 
of scalings will give rise to amplitude equations featur- 
ing the driving term of the same order of magnitude as 
the linear and nonlinear damping terms (so that the sta- 
tionary wobbling regimes become possible). We assume 
that the kink moves with a slowly varying, small velocity: 
v = tV where V = V(T l7 T 2 , . . .) is of order 1. 

Before embarking on the perturbation expansions, we 
transform Eq.JT]) to the co-moving reference frame: 

^(1 + p) 2 ^ TT - V(l + p)cb iT ~ y (1 + P)^ - 

—<j) + 3 = hcos(uioT)(p + yv<j>£ — 7(1 + p)<j> T - 

(3) 

Here 

£ = x- ( v(t')dt'. (4) 
Jo 

We have also changed t — ► r, where 
Qt = loqt. 
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We now expand the field 0(£, r) about the kink 0o = 
tanh£: 

= O + £01 + e 2 02 + ■ ■ ■ ■ (5) 

We also define "slow" space and time variables 

X n = e n £, T n = e n r, n 0.1.2 

wit h the standard short-hand notation d n = d/dX n , 
D n = d/dT n . Substituting {5| into the 4 equation 
([3]), making use of the chain-rule expansions d/d^ = 
d + edi + e 2 d 2 + . . ., d/dr = D + eD 1 + e 2 D 2 + . . ., 
and equating coefficients of like powers of e, we obtain 
a sequence of linear partial differential equations — jus t 
as we have done in the case of the free wobbling [2lj . 
As in the case of the undamped undriven 4 equation, 
the first-order perturbation is chosen to include only the 
wobbling mode, 

01 = AiXi, ...;T 1: ...)sechX Q tSinhX Q e luJoTo + c.c, (6) 

while the quadratic correction satisfies the partial differ- 
ential equation 

\Dl4> 2 +C^=F 2 {X ,...-T ,...) (7) 

with 

C = -\dl - 1 + 30o = -\dl + 2-3 sech 2 X 

and 

F 2 = (dodx - D D 1 )<j> 1 - 30o0? + VDadofa 

+^D 1 Vd <t> ~ \v 2 dU . (8) 

As in [21[ , the second-order perturbation is taken to con- 
sist just of the harmonics present in the forcing function 
©: 

2 = + + c.c. + <pWe 2i "° T ° + c.c. (9) 

Here the coefficient functions ip 2 (Xo), ip 2 (Xo), and 

(2) 

if 2 (Ao) are found by solving the corresponding ordi- 
nary differential equations. The solvability condition for 



F 3 = (dodi - DoD^ + (d d 2 - D Q D 2 )<i>i + \{d 2 
+^^001 + \D 2 Vd a cj> - \V 2 d 2 ^ - 



The cubic correction 03 consists of harmonics present in 
the function F3. The solvability condition for the first 



the first of these equations is D{V = and the solution 
is 

</4 0) = 2 2 sech 2 A tanh A 
+ ("T~ 3|A|2 ) X osech 2 A (10) 
(see [HI). The solution of the last equation is 

^ 2) = A 2 f 1 (X ), (11) 

with 

/i(A ) = I {6 tanh A sech 2 A 
+ (3 - tanh 2 A + ik Q tanh A )[J 2 * (A ) - J|°]e ife ° x ° 
+ (3 - tanh 2 A - ifcotanhAo)J 2 (Ao)e^ feQXo }. (12) 
Here the function J 2 (Xo) is defined by the integral 

Jn(X )= [ ° e ifco « sech" £ d£ (k = V8), (13) 
J — 00 

with n = 2. The constant J^ is the asymptotic value of 
J 2 (A ) as A — > 00: 

J~= lim J„(A ). (14) 

A|] — >00 

Finally, the nonsecularity condition associated with 
the equation for the coefficient function tp^ (Xo) is 
D\A — 0; with this condition in place, the solution 
ip 2 (Xo) is bounded for all Ao and decays as |Ao| — > 00. 
However, this decay is not fast enough [2l|; hence the 
term ip 2 1 \Xo)e l ' JJoTo has a gwasisecular behaviour at the 
infinities and has to be set to zero. This is achieved by 
imposing the condition [2l[ 

diA + iu VA = 0. (15) 

Proceeding to the order e 3 , we find the PDE 

§^03 + ^03 = ^3, (16) 

where 



- £ 2 )0! - 0? - 60o0i02 + VD d Q <f> 2 + VDodifa 

- rD o 0! + ry5 o 0o - RD 2 ^ + ±He^ oTo <p + c.c. (17) 

I 

harmonic gives the amplitude equation 

UhA + ^f\A\ 2 A -u (R + ^)A=^H- iTA, 

(18) 
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while the solvability condition for the zeroth harmonic 
produces 

D 2 V = -2TV. 

Letting a = eA and keeping in mind that D±A = and 
D\V = 0, we can rewrite these two equations in terms of 
the unsealed variables. This gives a system of two master 
equations 



-7a — iuopa + 



1 SI |2 1 2 



-1 — UjQtl 



0(\a\») 
2 1V + 0(\a\ 5 ) 



(19a) 
(19b) 



where the overdots indicate differentiation with respect 
to t and the complex coefficient £ was evaluated in [2l| : 



C = Cr + Ki 



-0.8509 + i 0.04636. 



(20) 



B. Reduced Two-Dimensional Dynamics 



Since a is complex, Eq. (|19p defines a dynamical system 
in three dimensions. However, equation (|19b[) will damp 



the variable v until it is of order 



and this will make 



the term v 2 a negligible in Eq. (|19a[ ). Thus, after an initial 
transient, the dynamics will be determined by the two- 
dimensional system (|19a[) with v = 0. Next, the natural 
wobbling amplitude a may depend, parametrically, on £. 
However, equation (|15[) and the fact that v — > as t 
grows, imply that a may only depend on £ via X2, X3, 
etc. That is, the dependence is weak. 
Letting a = re~ ie , Eq. 1(155]) yields 



r = — 2^ +7( r osm( 



r), 



(21) 



where ro = ^wo^/7- For all a with |a| > ro, the right- 
hand side of f2T]) is negative and so no trajectories can 
escape to infinity. On the other hand, applying Dulac's 
criterion (with Dulac's function equal to a constant), one 
can easily ascertain that Eq. (|19a[) with v = does not 
have closed orbits. Hence, all trajectories must flow to- 
wards fixed points with finite \a\. 

The fixed points of the system (jT9aJ) are given by the 
equation 



, , . v i^oC, ,2 nivo, 
(7 + iuj p)a — \a\ a = —h. 



From Eq. ((22]) . the absolute value of a satisfies 

H(\a\ 2 ) =h, 
where the function 7i(|a| 2 ) is defined by 



(22) 



(23) 



7T Z 



J_ 



7| a| 2 



Cr 1 r 



(24) 



Assume, first, that p > po, where 

1 O - V3Ca. 



pod) 



w Cr + y/K 



-7 



-I.I397. 



(25) 



In this case 7i(|a| 2 ) is a monotonically growing function, 
with the range (0, 00). Eq. (f23]) has a single positive root 
I a. 1 2 for any h and the dynamical system ((15]) has a single 
stationary point. This fixed point is always stable. 

Now let p < p . Here, the range of the function H(\a\ 2 ) 
is still (0,oo); however the function grows for small and 
large values of |a| 2 but decreases in the intermediate in- 
terval |a_| 2 < |a| 2 < |a + | 2 , where 



l2 2 2a ± J a 2 - 3(3 2 
|a±| 2 =3 ^ -> ( 26a ) 



CrP-Ci 



J_. 

UJQ ' 



(3 = (iP + (r 



J_ 

UJQ 



(26b) 



Consequently, Eq. pS]) has a single root for small and 
large values of h and three roots in the intermediate re- 
gion defined by h+ < h < h-, where 

h±=H(\a ± \ 2 ). (27) 

For the dynamical system ((T9)) this implies that there is 
only one fixed point (which is stable) for small and large 
h, but as h approaches the value h+ from below or the 
value h- from above, two new fixed points are born in 
a saddle-node bifurcation. The region h + < h < h- is 
characterised by bistability; an adiabatic variation of h 
will result in hysteretic transitions between two stable 
fixed points. 

The existence of hysteresis has been verified in the 
direct numerical simulations of Eq.([T]), with 7 = 0.01. 
Starting with h = we increased h past h—, and then 
reduced it back to zero. At each h-step, we used the final 
values of 4>(x) and <f>t{x) from the previous-step simu- 
lation as initial conditions for the new run. For each 
h we measured the value of a to which the numerical 
solution settled after transients died out. The result- 
ing amplitude \a\ is shown in FigfT] clearly visible is the 
hysteresis loop. FigfT] corresponds to simulations with 
p = —0.03; for p = —0.02 the hysteresis loop is smaller 
and for p — —0.01 it disappears completely. This is con- 
sistent with the value of po given by Eq. ((2"5"]) . The value 
h— = 0.008 at which the amplitude was recorded to jump 
from the bottom to the top branch in Fig[TJ and the value 
h + = 0.005 at which it dropped back as h was decreased, 
are also in agreement with the corresponding predictions 
of the amplitude equation. Namely, Eqs. p5]) . (IM]) and 
(|27]l give h- = 8.2 x 10~ 3 and h+ = 4.9 x 10" 3 . 

For any given h, equation (]23[) can be regarded as 
a quadratic equation for the detuning p where the co- 
efficients are explicit functions of \a\ 2 . There are two 



roots pi,2 for I a. 1 2 smaller than 



\a\ 2 > |a rcs | 2 , where |a r , 
the equation 



and none for 



is a unique positive root of 



64, 



_7_ 

U) 



si I |2 
— Ores 



= K 2 



(28) 



5 



0.35 




FIG. 1: The hysteresis loop observed in the 1 : 1 para- 
metrically driven (f> 4 equation, Eq.JTJ, with 7 = 0.01 and 
p = —0.03. The driving strength h is increased from to 
0.0125 in increments of 5 x 10 -4 , and then reduced back to 
as indicated by arrows. Crosses mark simulations of the 
equation fl}; continuous and dashed lines depict stable and 
unstable fixed points of the amplitude equation (|19a[l with 
v = 0. 

The value |a ros | defined by (|28j) gives the largest am- 
plitude of the wobbling achievable for the given driving 
strength h. The corresponding value of the detuning, 

Prcs = Ykes| 2 <0, (29) 

ensures the strongest resonance. As was expected, the 
strongest resonance is achieved with negative detuning. 

C. 1:1 Parametrically Driven Wobbler 

For large times, the asymptotic expansion for the 
damped-driven wobbler is given by Eq.(48) in [2l| where 
we just need to set v = and replace loq with Q: 

<t>{x,t) = tanh [(1 - 3|a| 2 )£] + a sech£ tanh£e if2t + c.c. 

+ 2\a\ 2 sech 2 f tanhf + a 2 /i(^)e 2t!! ' + c.c. + O (\a\ 3 ) . 

(30) 

Here £ = x — xq, where xq is a constant determined by 
initial conditions, and a is a stable fixed point of the dy- 
namical system (|19a|) with v — [a unique fixed point or 
one of the two stable fixed points depending on whether 
h is outside or inside the bistability interval (h+,h—)]. 
The function /i(£) is given by Eq. (fT2"|) . The interpreta- 
tion of different terms in (1301) i s the same as in the case 
of the freely wobbling kink [21] . 

Like the corresponding formula for the free wobbler, 
the expansion ([3H)) is only valid at distances |£| = 0(1). 
For larger distances one has to use the outer expansions 

= ±l + e 2 2 +e 3 3 + ..., (31) 



with coefficients <f> n determined as in section V of [2l|. 
The analysis of the outer equations produces results 
equivalent to those in [2l[: The second- harmonic radi- 
ation propagates away from the core of the kink at the 
group velocity, leaving in its wake a sinusoidal wave with 
the frequency 20, wavenumber fco = \/8 and constant 
amplitude of the order |a| 2 . 

Unlike the case of the free wobbling of the kink, the 
frequency of the oscillation is not determined by its am- 
plitude but is locked to the frequency of the driver, O. 
Another difference from the undamped undriven case is 
that the driven oscillations of the wobbler do not die out 
as t — > 00. Instead, the amplitude of the oscillations ap- 
proaches a nonzero constant value which is determined 
by the parameters of the damping and driving and - 
in the bistable region — by the initial conditions. On 
the other hand, the asymptotic velocity of the damped- 
driven wobbler is zero. 

It is interesting to note that unlike in the case of the 
parametrically driven damped linear oscillator [l8| or 
damped-driven breather of the sine-Gordon or </> 4 equa- 
tion [19| , there is no threshold driving strength in the case 
of the damped-driven wobbler. No matter how small is 
h, the amplitude a(t) will not decay to zero as t — > 00. 

III. 2 : 1 PARAMETRIC RESONANCE 

A. Asymptotic Expansion 

It is a textbook fact that the strongest parametric res- 
onance is achieved when the parameter of the oscillator 
is varied at double its natural frequency. With an eye 
to the detection of the most efficient driving regime for 
the wobbling kink, we now consider the driving frequency 
close to twice its natural wobbling frequency: 

\<i>tt 7& - [1 + h cos(20t)] <j> + <j> 3 = 0. (32) 

As before, 

n = uj {l + p), P = e 2 R, 7 = e 2 r, 

but now we use a different scaling for h: 

h = e 2 H. 

We transform the equation in exactly the same way as 
we did in the previous section; this yields 

i(l + p) 2 4> TT - w(l + p)<t>zr ~ y (1 + p)h - ^Y^^i 
— (j) + (j) 3 — h cos(2w o t)0 + 7«<?^ — 7(1 + p)(j>T- 

The perturbation expansion is unchanged from the un- 
damped undriven case at O (e 1 ). With the addition of 
the e 2 -strong driving, the equation at O (e 2 ) acquires ad- 
ditional terms on the right hand side as compared to 
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Eq.©: 

\Dl<f> 2 + C<i> 2 =F 2 {X ,...;T , 



Here F 2 is as in Eq.©. The zeroth- and first-harmonic 
components of <j) 2 are not affected by this extra term. 
Namely, assuming that the solution is in the form ([9]) 
and setting D{V — 0, we get Eq. fTO]) for ip 2 °^ while im- 
posing D\A = 0, Eq. (fT5|) produces (p 2 = 0. As for the 

(2) 

coefficient function ip 2 , we obtain 



Y2 



(2) 



A 2 f 1 (X ) + Hf 2 (X ), 



where the function f 2 (Xo) satisfies 



(C-6)f 2 (X ) 



tanh X, 



o- 



(33) 



(34) 



We note that the value 6 lies in the continuous spectrum 
of the operator £, and so in order to determine f 2 (X ) 
uniquely, one has to impose two additional conditions 
fixing the coefficients of two bounded homogeneous solu- 
tions that can be added to f 2 . We do this by requiring the 
absence of incoming radiation. The particular solution of 
([34]) which obeys these radiation boundary conditions is 



h(X a ) = -^/ 1 (X ) + ^tanhX (2scch 2 X -3), (35) 

where the function fi(Xo) is given by Eq. (fT2|) . In what 
follows, we will use the fact that J 2 (Xq) is an odd func- 
tion. 

The first term in the right-hand side of (|33)) describes 
the familiar second-harmonic radiation from the freely 
wobbling kink. The second term consists of the induced 
second-harmonic radiation and a standing wave — also 
excited by the forcing. 

At the order e 3 , we get the equation (flU)) where F 3 
is given by Eq. (fT7|) with the term 2 -He lu)oT ° 4>o replaced 
with ^He 2luJoTo (f>i. The amplitude equation for A, which 
arises as the solvability condition for the first harmonic, 
is now 

iD 2 A + ^f\A\ 2 A -u Q (^+R^A = ^HA* - iTA, 

(36) 



where 



/CO 
[77 sech 2 Xq tanh 2 X 
-00 

6 sech 2 X a tanh 3 X f 2 (X )] dX . 



The imaginary part of this integral is 
a/ = — = 0.003863; 

for the real part we find, numerically, 
a R = 0.5958. 



The solvability condition for the 0th harmonic gives 

D 2 V = -2TV. 

We finally write the amplitude and velocity equations 
in terms of the natural variables a = eA and v = eV and 
the unsealed time t (as in the previous section): 



a = —7a — iwopa + -ioJo(\a\ 2 a 

— -iujQV 2 a — -iujQaha* + O (|a| 5 ) 
v = -2 1V + O (\a\ b ) . 



(37a) 
(37b) 



B. Reduced Dynamics in Two Dimensions 

As in the previous case of the 1 : 1 parametrically 
driven wobbler, the velocity tends to zero as t — > 00 and 
the evolution of a(t) is governed by the dynamical sys- 
tem (|37aj) with v = 0. This two-dimensional dynamical 
system does not have periodic orbits, as one can read- 
ily check using Dulac's criterion. Letting a = re and 
a = |cr|e lArg<T , Eq.([37al) yields 

r = - 7 r + y [r 2 sin(20 + Arg a) - r 2 ] , (38) 

where r 1 = (\a\/(i)h. Since the right-hand side of 
is negative for all a with \a\ > ro, no trajectories can 
escape to infinity. Therefore, all trajectories should flow 
to one of the fixed points. The fixed points are given by 
the equation 

7a + lojQpa — 1— — \a\ a = —i—^-ha . (39) 

One fixed point is trivial, a = 0; this fixed point is stable 
if h < where 



1£L 
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hi = 1 7 1 - 2 



P 



(40) 



and unstable otherwise. For the nontrivial points, we get 
H(\a\ 2 )=h, (41) 



where 



H z = 



5/ I |2 



w 



Assume, first, that p > po, where 

po(l) = — 7^7 = -0.03146 7 . 



(42) 



The function 7i(|a| 2 ) with p in this parameter range is 
monotonically growing, from h + to infinity. The equa- 
tion (T4"T]) has one root provided h > h + , and no roots 
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otherwise. Consequently, in the region h > h+ the dy- 
namical system (|37p has 2 stable fixed points a\ and —0,%, 
where 



2a+ VlC^pfe 2 ~4/3 2 
ICI 2 



(43) 



with a and f3 as in (|26bp . In the region h < h + , the 
(stable) fixed point at the origin is the only fixed point 
available in the system. (Thus we have a supercritical 
pitchfork bifurcation as h is increased through h = h + .) 

Assume now p < po- As |a| 2 grows from zero to infinity, 
the function Ti decreases from h + to its lowest value of 
h-, where 



(44) 



and then increases to infinity. Therefore, for p < po ; the 
dynamical system (|37|) has one fixed point at the origin 
for small driving strengths < h < h-\ 5 fixed points 
ai,a2,0, — ai, and — a% for the intermediate strengths 
h- < h < h+; and 3 fixed points a l7 0, and — a\ for 
h > h + . Here h+ is given by Eq. (j4*D|) and /i_ by (|3S)). 
The nontrivial fixed points a\ and a 2 are born in a saddle- 
node bifurcation at h = h—. At h = a subcritical 
pitchfork bifurcation occurs; here, the point a 2 merges 
with the trivial fixed point. Therefore, out of the two 
nontrivial fixed points a\ and 02, the stable one is a%, i.e. 
the fixed point with the larger absolute value — given by 
Eq. (|4"3")) . In summary, for h < h- all trajectories flow to 
the origin; for h > h + they are attracted to the nontrivial 
fixed points ±ai, and, finally, in the region h- < h < h + 
we have a tristability between a = and a — ±ai. 

These predictions of the amplitude equation were 
compared to results of direct numerical simulations of 
Eq.(|32p, with 7 = 0.005. As in the previous section, 
we increased h past h+ and then reduced it to values 
under h-. Figf2] shows the hysteresis loop arising for 
p = —0.005. For p = —0.003 the hysteresis was less 
pronounced and for p = it was seen to disappear com- 
pletely. These observations are consistent with the value 
of the critical detuning (|42[) which, for 7 = 0.005, equals 
po = —1.573 x 10 -4 . The bifurcation values h± ob- 
served in simulations with p = —0.005 (h+ = 0.019 and 
h- = 0.011), are also in agreement with the predictions 
of the amplitude equation (which gives h + = 0.01938 and 
h- = 0.01059). 



C. 2:1 Parametrically Driven Wobbler 

Restricting ourselves to the t — > 00 asymptotic be- 
havior of <p(x,t), the leading orders of the perturbation 
expansion in the case of the subharmonic response are 

4>(x,t) = tanh [(1 - 3|a| 2 )C] + ascch^tanhfe in * + c.c. 
+ 2|a| 2 sech 2 ^tanh^ 



[a 2 fi(0 + hf 2 (0] 



2iQt 



+ C.C. 



0(\c 



0.16 



0.12 



0.08 



0.04 



0.008 



0.012 



0.016 



0.02 



0.024 



FIG. 2: The hysteresis loop in the 2 : 1 parametrically driven 
(j> 4 equation, Eq.fl32J, with 7 = 0.005 and p = -0.005. The 
driving strength is increased from 0.005 to 0.025 in increments 
of 8 x 10~ 4 , and then reduced back to 0.005. Crosses mark 
results of simulations of Eq. (|32[l . The continuous and dashed 
lines show the stable and unstable fixed points of the ampli- 
tude equation (|37ap with v = 0. 



Here £ = x — xq\ a is a stable fixed point (zero or nonzero) 
given by one of the roots of (13"9")) , and the functions /1 and 
ji are defined by Eqs.fH]) and (|35|) . respectively. The 
main difference from the case of the 1 : 1 parametric res- 
onance is that the amplitude of the wobbling approaches 
a nonzero value only if the driver's strength exceeds a 
certain threshold; this threshold value is given by h + in 
the region p > po and by h- in the region p < pq. If h 
lies below the threshold, the wobbling dies out and we 
need to set a = in Eq. (|45|) . We also note that the 2 : 1 
resonant driving excites a standing wave and radiation 
with the frequency 2Q and amplitude proportional to h 
[the /2-term in (TiS]) ]. 

As Eq. (|3"0)) of the previous section, the expansion (|4"5"|) 
is only valid on the lengthscale £ = O (1). To describe the 
waveform at longer distances, we need to invoke the outer 
expansions (|3"Tj) . Evaluating the coefficient of the term 
e 2 in these expansions and matching it to the "inner" 
expression Q33p in the overlap region, we obtain 



b 2 = ±JB±e 



i(2uj T TkoX ) 



+ C.C. 



H 



cos(2w T ), (46) 



(45) 



where the top and bottom sign pertain to the regions 
A > and X < 0, respectively. In Eq.lgBJ), J = (2 - 
ik )J?> and the functions B± = B±{X X , X 2 , T x , T 2 , ...) 
satisfy 

B ± (0,0,...;Ti,T 2 ,...) = A 2 (0,0,...;T2,T 3 ,...)-g. (47) 

Eq. (|47p represents the boundary condition for the am- 
plitudes B±: the equations of motion for these variables 
arise at the order e 3 and coincide with Eqs.(33) of [2l| . 
The solution of these equations with the boundary condi- 
tion (|47p is qualitatively similar to the solution with the 
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"undriven" boundary condition B±(Q, 0, Ti, T 2 , ...) = 
A 2 (0, 0, T 2 , T 3 , ...). Namely, we have two outward- 
propagating waves leaving the amplitudes B± equal to 
the constant A 2 — H/96 in their wake. 



IV. HARMONIC VS SUBHARMONIC 
PARAMETRIC RESONANCE: QUALITATIVE 
COMPARISON 



With the amount of detail that we had to provide 
to justify our conclusions and derivations, the resonance 
mechanisms of the driven wobbling kink may not be easy 
to crystallise. The purpose of this short section is to dis- 
cuss the two parametric resonances qualitatively — in 
particular, to comment on their atypical hierarchy. 

We observed that the amplitude \a\ is of the order /i 1 / 3 
in the case of the harmonic resonance (i.e. the resonance 
arising when the driving frequency oj^ is near lo ) but 
only O (h 1 / 2 ) in the case of the subharmonic resonance 
(the resonance arising for LOd ~ 2loq). Thus the har- 
monic resonance is stronger than the subharmonic one, 
and this is precisely the opposite behaviour to what we 
might naively expect based on our intuition about the 
parametric driving. 

To explain this surprising behaviour qualitatively, we 
write the term hcos(nflt)(f> in Eqs.([T]) and (|32[) as 
hcos(nttt)(j)Q plus terms of order ha and smaller. This 
representation reveals that what was introduced as a 
parametric driver is, to the leading order, an external 
(direct) driving force. This driving force is nonhomogc- 
neous, i.e. its magnitude and direction vary with the 
distance, and it has odd spatial parity. When n = 1, the 
frequency of this driving force coincides with the natural 
frequency of the wobbler and its spatial parity coincides 
with the parity of the wobbling mode (which is also odd). 
As a result, we have a strong direct resonance. 

When n — 2, the external force h cos(2fii)0o is not 
in resonance with the wobbling frequency However the 
function hcos(2flt) acts as a parametric driver on the 
next term in the expansion of </> — that is, the product 
h cos(2Sli)e</>i has the resonant frequency. Importantly, 
this term has the "correct" , odd, parity as a function of 

In addition, the odd-parity force h cos(2flt)4>o gener- 
ates odd-parity radiation and the odd-parity standing 
wave, both with the frequency 20. This radiation and 
standing wave also couple to the wobbling mode, via the 
term e 3 



in Eg. (f52"| . This constitues a concurrent 
driving mechanism. Since each of the two mechanisms is 
indirect (i.e. requires the wobbling mode as a mediator 
for the frequency halving) and since the resulting effec- 
tive driving strength is proportional to the amplitude of 
the wobbling mode (assumed small), the response to the 
frequency 2f2 is weaker than to Q. 



V. 1:2 DIRECT RESONANCE 



A. Perturbation Expansion 

We start with the direct driving at half the natu- 
ral wobbling frequency. This and the following case of 
the 1 : 1 direct resonance were previously considered by 
Quintero, Sanchez and Mertens @, H[ and so we will be 
able to compare our results to theirs. The equation is 

\<t>tt - \<t>xx + l<t>t - <t> + 3 = h cos , (48) 

where, as in the previous sections, £1 = u>o(l + p). As 
before, we change the time variable so that fit = luqt 
and transform the equation to the moving frame: 

|(i + P ) 2 4> TT - w(i + p)^ T - y (l + p)4>i - ^Y^ka 

-cj) + <j> 3 = hcos(^j T ) + «70e- 7(1 + P)^r- (49) 

Keeping our standard scalings for the small parameters 
7 and p, 

1 = e 2 r, p = e 2 R, 
we choose a fractional-power scaling law for h: 
h = e z ' 2 H. 

This scaling will be shown to produce a balance of damp- 
ing and driving terms at the leading order in the ampli- 
tude equation. Expanding <f> in powers of e 1 / 2 , 

= 0o + e<£i + e 3/2 3 /2 + e 2 02 + e 5/2 ^ 5 /2 + • • • , 

where (f>o = tanhXo, and substituting in (|49[) . we obtain 
Eq. ([6]) for 4>\ . The partial differential equation arising at 
O U z ' 2 ) is 



-D, 



0^3/2 +£03/ 2 



H 



e ^ /2)T + cc 



The solution (/> 3 / 2 to this equation has the form 3 / 2 = 



v 3 3 /2 2 ' el< '" ^ 2 ' To + c - c -' wnere the coefficient func +; "-" •^ 1 / 2 " > 
satisfies the linear nonhomogeneous equation 



'-Ik? 1 -!- 



Since | is not an eigenvalue of the operator £, this equa- 
tion has a unique bounded solution. To determine it, we 
note that two homogeneous solutions of this equation, i.e. 
solutions of (£ — 3/8)y = 0, are given by Segur's formula 



1 



JpX Q 



x (2 - p 2 - 3ip t&nhXo - 3sech 2 X ) (50) 
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with p — ±i-\/l3/4 20]. Using these in the variation of 
parameters, we obtain 

^3/2 2) = T3 H ( l - 8 sech2 X o)- 

The term e l ^°/ 2 ^ To in the expansion of the wob- 

bling kink represents the background stationary wave in- 
duced by the driver. 

The equations arising at O (e 2 ) are the same as for the 
free wobbler and the 1 : 1 parametric resonance, Eqs.([7])- 
([5]). Hence the coefficients of the harmonic components 
of cj>2 are the same as in the undamped, undriven case. 
Namely, imposing the solvability conditions D{V = 
and DiA = 0, we obtain fHfy for <^ 0) and (HJ) for ipf>. 
We also impose Eq.(TJ5]) to obtain ^4 = 0- 

At the order e 5 / 2 we have the equation 

h D o4>5/2 + £<t>5/2 = -60O'/ , 103/2 + U£>o<9 </>3/2- 

Its solution consists of the ^th and |th harmonics with 
the coefficient functions 

= HAu a (X ) + iu HVu b (X ), 

tpffi = HAu c (X ), 

respectively. Here the functions u a , u b and u c satisfy 



and 



u a (X ) = -#(1 -8sech^A )sechA tanh^ A , 



27> 

8 / 



Uc(Xo) 



32 



sech AotanhAo, 



(1-8 sech 2 X ) sechXo tanh 2 X . 



In order to determine u c {Xq) uniquely, we impose the 
radiation boundary conditions. (These are necessary be- 
cause the value ^ lies in the continuous spectrum of the 
operator C.) The functions u a and u c are even, while u b 
is odd. These three functions can be easily found by solv- 
ing the above nonhomogeneous boundary- value problems 
numerically. 

Proceeding to the order e 3 , we find the equation (fl"6l) . 
where F 3 is given by Eq.(fTT|) with the term ^He luJ ° To ^Q + 
c.c. replaced with — 3(f>Q(j> 2 ^ 2 . The solvability conditions 
for this equation are 



D 2 V = -2TV 



(51) 



for the 0th harmonic, and 



D 2 A = -TA-iu RA+i(uo\A\ 2 A~^ujoV 2 A+-^ § iujoTrH 2 

(52) 

for the first harmonic. The latter equation includes both 
the damping and driving terms and so the resulting mas- 
ter equations could be expected to capture the essentials 
of the nearly-stationary wobbling of the kink (i.e. wob- 
bling in the vicinity of the fixed point of the amplitude 



equations which arises due the balance of the damping 
and driving terms). However the description provided by 
these amplitude equations — while being qualitatively 
correct — turns out to be insufficiently accurate when 
compared to numerical simulations of the full partial- 
differential equation (|4"8")) . (The source of this inaccuracy 
will be clarified below.) In search of greater accuracy, we 
shall proceed to higher orders. 

The solution of Eq. (TTHJ) has the form 

+ 4 2 >e 2iw ° T ° + c.c. + ^ 3) e 3MoTo + c.c. 

(2) 

The function tp^ ' is calculated to be zero; the coefficient 
of the 0th harmonic is given by 



(0) _ 16 tt2 



ji|H 2 (45Xo sech^ X - 3 tanhX 



- 128 sech 2 X tanh X ), 
and the one for the first harmonic component is 



(1) - -d 2 AX sech X tanh X + \A\ 2 Au d (X Q ) 



\iLU {T + iu R) (3 - A sech 2 X a ) + H 2 u e (X a ), (53) 



where the functions Ud{Xo) and u e (Xo) are the bounded 
solutions of the following nonhomogeneous equations: 

(C — ^)ud = |C sech Xq tanh Xq + 6 sech Xq tanh 2 Xq 



3X sech 2 X - f sech 2 X tanh X - f t (X ) 



(£- §H = - -^H 2 tanhX (l - 8scch 2 X ) 2 

180, 



(54) 



+ 7r sech Xq tanh Xq . 

Since (£ — 3/2) is a parity-preserving operator while 
the right-hand sides of the above equations are given 
by odd functions, and since the homogeneous solution 
y w = sech X tanh X is also an odd function, it follows 
that the nonhomogeneous solutions Ud and u e are both 
odd. This is the only fact about Ud and u e that we will 
need in this section — we do not need to know any detail 
of these functions here. Nevertheless, we do evaluate the 
solution Ud as it will be required later on, in the study 
of the 1 : 1 directly driven kink (section [VlJ) ; we evaluate 
it using the variation of parameters and numerical inte- 
gration. The nonhomogeneous solution is defined up to 
the addition of an arbitrary multiple of y w ; however this 
extra degree of freedom is fictitious as it can always be 
eliminated by a suitable rescaling of e. [Accordingly, the 
extra term proportional to y w cancels in the integral r\ 
where it appears in section IVII and does not contribute 
to the amplitude equations (|75p]. 

To eliminate the quasisecular term proportional to 
Xq sech Xq tanh Xq in (fS"3"]) , we set d 2 A = 0. 

It will not be necessary to calculate the third harmonic 

(3) 

component, <p\ , as this does not contribute to the 0th 
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or 1st harmonics at O (e 4 ) , and hence does not affect 
the e 4 -correction to the amplitude equations. Similarly, 
we shall not calculate (f>7/ 2 as it only contains fractional 
harmonic components which cannot impinge on the am- 
plitude equations at O (e 4 ) . Hence we skip the order 
e 7 / 2 . 



At O (e 4 ) , we obtain 

\D 2 Q( j> i + L^ = F A , (55) 

where 



F 4 = (fibft - D Q D{)<fo + (d Q d 2 - D Q D 2 )fo + - Dl)fo + (d Q d 3 - D Q D 3 )fa + {did 2 ~ DiD 2 )^i 
-3</> 2 </> 2 - 6^,01^3 - 30o^2 + VDodofo + VDodxfc + VD d 2 fa + VD 1 d 4> 2 + VD 1 d 1 <t> 1 + VD 2 d^ x 
+\lhVd Q (h. + ~D 3 Vdo<t> - \v 2 dl4> 2 - V 2 d d 1 cf> 1 - TD^ 2 - YD\4>x + TVdofa - RD 2 fo + VRD d fa. 

I 



The corresponding solvability conditions are 

D 3 V = (56) 

and 

D 3 A = -\iuj Q \H 2 A (57) 

where 

/oo 
sech X tanh X a [- ^ (45X sech 2 X 
-00 

— 3 tanh Xa — 128 sech 2 Xq tanh Xq) sech Xa tanh 2 Xq 

- f| tanhA (l - 8 sech 2 X )u a (X ) 

- f| tanhA (l - 8sech 2 X )w c (A ) 

- ^ sech X 1 anh X ( 1 - 8 sech 2 X ) 2 ] dX n . 

Numerically, 

X = Xr + i\j = -7.4656 - il.6785. 

Expanding the derivative d/dr as Dq + eD\ + e 2 D 2 + . . . 
and recalling that dr/dt = 1 + p, we combine equations 
([52)) and (JSTI). We also combine (J5TJ> and (J56]). This 
yields a system of two master equations: 

. - w oCi |2 2 

a = —7a — iLOopa + i— — \a\ a — % ~^ v a 

+ i-^TTLU h 2 - ±iLu Xh 2 a + O(\a\ 5 ) , (58a) 
v = -2jv + (|a| 5 ) . (58b) 

It is essential to combine the slow-scale equations in this 
way, rather than solving the individual equations with 
the assumption that the different scales are independent. 
Solving individual equations separately would be illegiti- 
mate because in integrating the equations one is covering 
more than one timescale. For example, solving Eqs. (|56p 
and (|57p we would be integrating over the scale T 3 which 
includes a shorter timescale T 2 . 




0.01 1 ■ ' ■ ' 1 

0.005 0.006 0.007 0.008 0.009 0.01 
h 

FIG. 3: The hysteresis loop in the 1:2 directly driven (j> 4 equa- 
tion, Eq.fl48}, with 7 = 0.001 and p = -0.002. The driv- 
ing strength h is increased from 0.005 to 0.01 in increments 
of 0.0002 and then reduced back to 0.005 (as indicated by 
the arrows). Crosses mark results of simulation of the PDE, 
Eq. (|48p . The continuous and dashed lines depict the stable 
and unstable fixed points of the amplitude equation (|58ap 
with v — 0. 



All terms in the right-hand side of Eq. (|58a[) are of the 
order |a| 3 , except the last term which is O (|a| 4 ). This 
last term is the correction coming from the fourth order 
of the perturbation expansion. As we have already men- 
tioned, the amplitude equations (|58[) without this term 
produce an inaccurate description of the dynamics in the 
region of interest (i.e. in the vicinity of the fixed points). 
On the other hand, if the above fourth-order term is in- 
cluded, the predictions of the amplitude equations (|58"|) 
turn out to be in good agreement with results of the di- 
rect numerical simulations of the full partial differential 
equation (|4"8")1 — see FigJ3] The substantial improvement 
in accuracy is due to the large value of the coefficient A. 
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B. Reduced Two-Dimensional System 



Introducing the notation 

/ Aj , o I Aft i 9 , I 480 , o 

7=7-y^o/j 2 , p' = P+^-h 2 , h' = - — h 2 , 
the amplitude equation (|58a[) can be written as 

/ • . I , -C 1 |2 * 2 

a = —7 a — zciJoP a + ^^l ! a — 2 UJ ° V a 

-i^u ti + O (M 5 ) , (59) 
o 

which has the same form as the amplitude equation for 
the 1 : 1 parametric resonance, Eq. (|19ap . Consequently, 
the dynamics of the 1 : 2 directly driven wobbling kink 
will have some similarities with the dynamics of the wob- 
bler driven by the 1 : 1 parametric force. 

According to Eq. (|58b[) , the velocity will be damped 
until it is so small [O (|a| 3 )] that it can be disregarded 
in Eq. (|58ap : hence after an initial transient the dynam- 
ics will be governed by Eq. (|58ap with v = 0. Simi- 
larly to Eq. lTOl) with v = 0, equation (|58ap does not 
have closed orbits. All trajectories crossing the circle 
\a\ = j^TTLdoh 2 / '7 flow inwards and so no trajectories 
can escape to infinity. Therefore, all trajectories must 
flow towards fixed points. If a is a fixed point, the abso- 
lute value la I satisfies 



4 
h 2 



|A| 2 



\a\ 2 I 169 



A M ,-i C BH 2 )-A,(-l + icz|a| 2 



P-\CR\a\ 2 



uj Z 



0. (60) 



The left-hand side of this equation is a bi-quadratic in 
h. Solving for h we obtain an explicit expression for 
h = ft-i,2(|a| 2 ); the roots \a\ 2 are found by inverting these 
explicit functions for each p, 7 and h. The lower branch 
of the bi-quadratic is plotted in figure [3] along with re- 
sults from the numerical simulations of the full PDE (|48f . 
We note that, as in the 1 : 1 parametrically driven <fi 4 
equation, there is no threshold driving strength for the 
existence of the nonzero wobbling amplitude here. 

Here it is appropriate to recall that p and 7 were as- 
sumed to be of the same order while h is O (7 3 ^ 4 ) . As 
in the case of the system (|19ap . the absence or presence 
of hysteresis in the dynamics depends on whether p is 
above or below the critical value ^0(7) given by Eq. (j25|) . 
If the difference p — po is positive and of order 7, there is 
only one root \a\ for each value of h. The corresponding 
fixed point is obviously stable. If, on the other hand, the 
difference p — po is negative (but still of the order 7), we 
have three roots \a\ 2 for each h — O (t 3 / 4 ) in the interval 
(h+,h-). These roots correspond to three fixed points, 
two of which are stable (see FigJ3]). The values h + and 



h- at which the subcritical bifurcations occur, are given, 
approximately, by 



h\ = I ( P 
Co 



2c Q d 



(61) 



where 



Q± = — 

V± 27 



P± = — [ci(3/3 2 - 5a 2 ) + 6a5 
18 

T(4cia-3(5) v / a 2 -3/3 2 



a (a 2 + 9/3 2 ) T (a 2 - 3/3 2 ) 3 / 2 



Co 



V 169 



ICI 4 = 0.3280, 



ci - (r Xr + C/Ai - 6.2747, 



<5-|C| 2 ( XrP-Xi — 

UJQ 



and a and (3 are defined by Eq. (|26bp . For 7 = 10 3 
and p = —2 x 10~ 3 , the bifurcation values obtained from 
Eq-dSH) are h + = 6.74 x 10" 3 and = 7.78 x 10" 3 
while the numerical simulations of the full PDE give 6.6 x 
10" 3 < h+ < 6.8 x 10~ 3 and 7.8 x 10~ 3 < h- < 8.0 x 
10~ 3 . For p = —3 x 10 -3 , the simulations show a more 
pronounced hysteresis loop whereas for p = —1 x 10~ 3 , 
the hysteresis was seen to disappear. (In both cases 7 
was kept at 10 -3 .) These observations are consistent 
with the value of po given by the amplitude equations. 
[For 7 = 10~ 3 , Eq.© gives p = -1.1 x 10~ 3 .] 



C. 1:2 Directly Driven Wobbler 

Finally, we produce the first several orders of the per- 
turbation expansion for the 1 : 2 directly driven wobbling 
kink: 

4>{x,t) = tanh[(l - 3|a| 2 )£] + a sech £ tanh £e i0t + c.c. 
+ *Lh(l - 8 sech 2 eyW 2 )* + c . c . 

+ 2|a| 2 sech 2 £ tanh £, + a 2 h{£)e 2mt + c.c. + O (\a\ b/2 ^ 



(62) 

When t is sufficiently large, the variable £ in this expres- 
sion equals x — xq (where xq is a constant determined 
by the initial conditions) and a is a stable fixed point of 
the system (|58ap with v — [a unique fixed point or one 
of the two stable fixed points depending on whether h is 
outside or inside the bistability interval (h + , h-).] 

The interpretation of terms in (|62p is the same as in the 
previous sections. The frequency of the wobbling [where 
the wobbling mode is given by the sum of the third and 
second term in the first line in l]62p] is locked to double 
the driving frequency. The term proportional to h in the 
second line describes a stationary wave induced by the 
driver. As in the previous sections, the expansion 
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only valid at the length scale |£| = O (1). The standard 
analysis involving outer expansions demonstrates that for 
larger distances, we have groups of second-harmonic ra- 
diation waves moving away from the kink and leaving in 
their wake a sinusoidal waveform of constant amplitude. 

D. Qualitative Analysis 

The driving term h cos(^t) is not in resonance with the 
natural frequency of the wobbler, nor does its parity coin- 
cide with the parity of the wobbling mode. Therefore the 
ability of the 1 : 2 direct driving to sustain the wobbling 
is surprising and requires a qualitative explanation. 

The authors of |2j propose that the mechanism which 
brings about the unexpected superharmonic resonance is 
the coupling of the translation mode and the wobbling 
mode. Our explanation for this phenomenon is rather 
different and unrelated to the translation mode. The 
way the driver affects the wobbler is by exciting an even- 
parity standing wave (^3/2) at the frequency Q/2 which 
then undergoes nonlinear frequency doubling and par- 
ity transmutation through the term e 3 (j)o(j) 2 , 2 in Eq. (|48[) . 
This latter term serves as an effective driver to the wob- 
bling mode; it has the resonant frequency and "correct" 
parity. 

Since this mechanism involves a two-stage process and 
the resulting effective driving strength is proportional to 
h , this type of driving produces a relatively weak re- 
sponse. 

E. Chaotic Wobblers? 

The authors of Ref.0, [|[ observed chaotic kink dy- 
namics in numerical simulations of the 1 : 2 directly 
driven wobbling kink. An indirect confirmation of 
the existence of chaotic motions comes also from the 
collective-coordinate approach which predicts an un- 
bounded growth of the kink's width, energy and velocity 
at resonance 0, On the other hand, our amplitude 
equations (|58[) with 7 7^ reduce to a two-dimensional 
dynamical system which can obviously not exhibit any 
chaotic attractors. 

To find an explanation for this disagreement, we have 
carried out a series of numerical simulations of the partial 
differential equation (|48[) at a range of driving strengths 
and damping coefficients. In all our experiments, we con- 
fined ourselves to zero detuning, p = 0. We could not 
detect any sign of chaotic dynamics for h smaller than 
a certain minimum value, not even in the undamped 
case. However for h greater or equal than 0.05 and suf- 
ficiently small 7, our numerical simulations did reveal 
kinks performing erratic motion, where initially close pro- 
files were seen to diverge exponentially fast. For h = 0.05, 
0.06 and 0.08, chaos was observed in simulations with 7 
smaller or equal to 10~ 3 , 2 x 10~ 3 , and 6 x 10 -3 , respec- 
tively, whereas the same sequence of h values paired with 



7 = 2 x 10~ 3 , 3 x 10~ 3 and 7 x 10" 3 , respectively, did 
not feature any chaotic trajectories. Therefore chaotic at- 
tractors may only arise when the damping is extremely 
weak, much weaker than O (/1 4 / 3 ) . This is the reason 
why the chaotic dynamics is not captured by our am- 
plitude equations (|58p which have been derived on the 
assumption that h — O (e 3//2 ) and 7 = O (e 2 ). 

The description of chaotic motions by means of ampli- 
tude equations is a topic of future research. We expect 
our asymptotic method to remain applicable in this situa- 
tion, with the appropriate adjustment of the scaling laws 
of the a and v variables and parameters of the damping 
and driving. 

VI. 1:1 DIRECT RESONANCE 
A. Multiscale Expansion 

Finally, we explore the effect of the direct driving near 
the natural wobbling frequency of the kink. The equation 
is 

\<\>u - \4>xx + l4>t - <t> + <j? = h cos(fit), (63) 

where 57 = u>q(1 + p). We let a = eA and adopt the 
following scalings for the three small parameters: 

h = eH, 7 = e 2 r, p = e 2 R. (64) 

This time, we assume that the velocity is scaled as 
v = e 2 V (and not as v — eV). We shall find a non- 
trivial evolution equation for v, and with the above scal- 
ings, the leading-order dynamics of a and v will occur 
on the same timescale. While other scalings could be in- 
vestigated, the variables v and a would then change on 
different timescales and so would effectively be decoupled 
for small e. Therefore, the chosen scalings correspond to 
the richest, three-dimensional, dynamics. Rescaling the 
time so that Qt — luqt and transforming to the moving 
frame as in Eq.(U]), the equation (|6"3")) becomes 

|(i + P ) 2 4> TT - w(i + p)hr - y (1 + p)fc - ^Y^ka 

- 4> + (j) 3 = h cos(w t) + V704 - 7(1 + p)4> T - (65) 

We expand as in Eq.lJS]). 

With the driving amplitude of the order e, the linear 
perturbation consists of the wobbling mode and a stand- 
ing wave excited by the driver: 

0i = [AsechXotanhXo + H(l - 2sech 2 X )] e lu ° T ° + c.c. 
Next, at O (e 2 ), we obtain the equation ([7]) with 

F 2 = (d d 1 - DqD^x - 3<M> 2 . 

Once transients have died out, the solution to Eq.([7]) 
will consist only of the harmonics present in the forc- 
ing, i.e. will have the form The solvability condition 
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for the first-harmonic component is D\A = 0; assum- 
ing that this condition is in place, we obtain tp^^Xo) = 
— diAXo sechXo tanhXo. To avoid the quasisecular be- 
haviour at infinity we impose d\A — 0, which results in 

<pW = 0. (66) 

The other two harmonic components of the quadratic 
correction fa have the coefficients 

ip 2 0) = \A\ 2 sech 2 X (2tanhX - 3X ) 
+ H 2 (9X a sech 2 X - 3 tanh X Q - 8 sech 2 X Q tanh X ) 
-4iJ(A + A*)sechX (l + sech 2 X ) (67) 

and 

<p {2) = A 2 A(X ) + AHf 3 (X Q ) + H 2 U(X ), (68) 

where f\ is as in Eq. (|12p , and the functions fa and f± are 
defined by 

f 3 (X ) = \ scchX - 4 sech 3 X 

- §ik (3 - tanh 2 X a + ik tanh X ) [J x * (X ) - J°°]e ifeo *° 

+ ||ife (3 - tanh 2 X - ifc tanhXo) Ji(X )e" !fcoJfo 

(69) 



and 



U(X ) = -\h(*o) + itanhX (3-2sech 2 X ). (70) 



The function Ji(Xq) is given by the integral (flU)) with 
n = 1. One can show that Js^Xq) is an even function 
and f<i{Xo) is odd. 



We note a quasisecular term (9H 2 — 3\A\ 2 )Xo sech 2 Xq 
in Eq. (|67p : this term does not lead to the nonuniformity 
of the expansion as it can be incorporated in the variable 
width of the kink. [Sec Eq.® below.] 



The partial differential equation arising at the order 
e 3 , is Eq. (fT6|) . with F 3 given by 



F 3 = (dodi - DoD^fa + (d d 2 - D D 2 )fa + ±(d 2 - D 2 )fa - 4>\ 
-Gfafafa + VD Q d fa + \D x Vd Q fa - TD fa - RD 2 fa. 

I 

The solvability conditions give rise to amplitude equa- for the first harmonic. In Eq. (|72p . we have introduced 
tions: 

£>iV = (71) 
for the zeroth harmonic, and 

D 2 A + TA + iuj RA - iiw C| A\ 2 A - \WH 

+ \iuj vH 2 A + ^iLu fiH 2 A* = (72) 



/OO 
sechXotanhXo [— 6tanhXo(l — 2 sech 2 Xo)f^(Xo) — 6 sech Xq tanh 2 Xq(9Xq sech 2 Xq — 3tanhXg 
-oo 

-8sech 2 X tanhX ) + 24 tanh X (sech X + sech 3 X )(l - 2 sech 2 X ) - 6sechX tanhX (l - 2 sech 2 Xq) 2 ] dX Q , 
and 

/oo 
sechX tanhX [24tanhX (l - 2sech 2 X Q )(sechX + sech 3 X ) 
-oo 

— 3 sechXo tanhXo(l — 2sech 2 Xo) 2 — 6sechXo tanh 2 -Xq/^Xo)] cLXq. 
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Numerically, 

v = 4.159 - i0.3258, fi = 1.022 + i0.1623. 

We note that the velocity enters the amplitude equa- 
tion (|72")l as a coefficient in front of one of its two driving 
terms. On the other hand, Eq. dTTI) implies that V does 
not tend to zero — at least on the timescale T\ . In order 
to check whether the velocity decays on a longer timescale 
and hence whether the translational motion can drive the 
wobbling, we take the expansion to higher orders. 

The cubic correction has the form 

3 = <4 0) +4 1) e iwoTo +c.c. 

+ 4V iWoT ° + c.c. + ^ 3 3) e 3 ^ To + c.c. 



(£-3/2)ui 
(£ - 3/2)u 2 

(£ - 3/2)u 3 
(£ - 3/2)u 4 
(£ - 3/2)u 5 



(£ - 3/2)u 6 

and 

(£-3/2)u 7 



Like the functions ua and u e of the previous section, the 
solutions u n (Xo) are defined up to the addition of a mul- 
tiple of y w . As in the previous section, this does not 
provide any extra degrees of freedom and the multiple of 
y w cancels out in the integrals r] and \ below. The solu- 
tions Ui, i*5, and ue are odd, while u 2 , U3, 114 and U7 can 
be chosen to be even functions. The only fact about u\ 
that we need is that it is a real solution; owing to its re- 
ality, u\ does not contribute to the solvability conditions 
below. The solutions u n (Xo) with n = 2, 7 are deter- 



The 0th harmonic component, , is evaluated to be 
zero, and the coefficient of the first-harmonic component 
is 

<p^ - \A\ 2 Au d (X ) 

- (d 2 A + iuj VA)X sechX tanhX 
+ iu) VH Ul (Xo) + H\A\ 2 u 2 (X ) + HA 2 u 3 (X ) 
+ HA 2 u 4 (X ) - |io;o(r + iuj R){3 - 4scch 2 X ) 
+ H 2 Au 5 (X Q ) + H 2 A*u 6 (X ) + H 3 u 7 (X Q ). (73) 

Here u^Xq) was defined in the previous section as the 
bounded solution of Eq. (|54"|) , and the functions u n (X ) 
(n = 1, ...,7) are the bounded solutions of the following 
nonhomogeneous equations: 



mined using the variation of parameters and numerical 
integration. 

To eliminate the quasisecular term proportional to 
X sech X tanh X in ([75)) . we set d 2 A = —iuj VA. 

It is not necessary to calculate the second and third 
harmonic components, ip^ and ip^ , as these do not con- 
tribute to the zeroth harmonic at fourth order in e, where 
the leading order behaviour of V will reveal itself. The 
equation arising at C(e 4 ) is Eq. ([55]) . with 



I 

37T 

= — 4sech Xq tanh Xq + — sech Xo tanh Xo, 

= -6 sech 2 X Q tanh 2 X (l - 2 sech 2 X ) + 24 sech 2 X tanh 2 X (l + sech 2 X ) 

— 6 sech 2 X Q tanhX (2 tanhX - 3X )(1 - 2 sech 2 X ) ~ 6 sechX tanh 2 Xo/ 3 (X ), 
= 3 sech 2 X tanh 2 X (7 + 10 sech 2 X ) + 12 sech 2 X tanhX /i(Xo), 

= -6tanhX / 1 (X ), 

= -6tanhX (l - 2 sech 2 X )/ 3 (X ) + 24 sech X tanh X (l + sech 2 X )(l - 2 sech 2 X ) 

— 6 sechXo tanh 2 Xo(9Xo sech 2 Xq — 3 tanhXo — 8 sech 2 Xo tanhXo) 

— 6 sech Xq tanh Xq (1 — 2 sech 2 Xq ) 2 — § v sech Xq tanh Xq , 

= -6sechX tanh 2 Xo/ 4 (Xo) + 24 sech X tanh X (l + sech 2 X )(l - 2 sech 2 X ) 

— 3sechXotanhXo(l — 2sech 2 Xo) 2 — |/x sech Xo tanhXo, 

= -6tanhX (l - 2 sech 2 X )/ 4 (X ) - 3(1 - 2 sech 2 X ) 3 

— 6 sech Xo tanh 2 Xo (9Xo sech 2 Xo — 3 tanh Xo — 8 sech 2 Xo tanh Xo ) . 

I 



I 

F 4 = (dodi - A)£i)03 + (dvd 2 - DoD 2 )fa + \{dl - D\)<p 2 + {d d 3 - £> o A0& + [d x d 2 - D x D 2 )<t>i 
-30 2 2 - 6000103 - 30o0^ + VD do<h + ra O 5i0i + ^D 2 Vd ^> - ~V 2 d$cf> - RD 2 fo - TDafa + IVc\,0 o . 
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The solvability condition for the zeroth harmonic yields 

D 2 V = -2TV - lrjH\A\ 2 A + c.c. - ^xH 3 A + c.c. + ^ioj H{T - icu R)A + c.c, 



where 



and 



/oo 
sech 2 X Q [24sech 2 X tanh 2 X (sechX Q + sech 3 X ) - 3sech 2 X Q tanh 2 X f 3 *(X ) 
-CO 

-6sech 3 X tanhX (2tanhXo - 3X )(1 - 2sech 2 X ) - 6 sechX tanh X (l - 2 sech 2 X )fi(X a ) 

—6 sechX tanh 2 X U3(X a ) — 6 sechX tanh 2 X ua(Xq) 
-6sechX tanh 2 Xow 2 *(^o) + 24 sech 2 X Q tanhA (2 tanhX - 3X ) 
x (sech X + sech 3 X ) - 6tanhA (l - 2 sech 2 X )u d (X a ) - 6 teaihX f 1 (X )f 3 *{X Q )] dX , 



X = 



/OO 
sech 2 X [-6sechXotanh 2 X u 7 *(Xo) - 6tanhX (l - 2 sech 2 X )u 5 (X ) 
-OO 

-6tanhX (l - 2 sech 2 X )u 6 *(X ) + 24(1 - 2 sech 2 X ) 2 (sechX + sech 3 X ) 
+ 18 sech X tanh X (l + 2 sech 2 X ) x (9X Q sech 2 X - 3tanhX - 8 sech 2 X tanhX ) 
-3(1 - 2sech 2 X ) 2 / 3 (X ) - 6tanhX / 3 pf )/4*(X ) - 6 sechX tanhX (l - 2 sech 2 X )/ 4 *(X )] dX . 



(74) 



Numerically, 

?? = -2.005 - i0.3823, x = -12.21 - i0.5706. 

Writing a = eA as before, and combining D\A = 
with Eq. ([72|) with the help of the chain rule, we obtain the 
amplitude equation in terms of the unsealed parameters: 



a = —7a — iujQpa + ^iujQ(\a\ 2 a + ^-vh 

— ^iujQi'h 2 a — ^iLUQfih 2 a* + O (|&| 5 ) 



(75a) 



Similarly, combining Eq. ([7T|) with ([71)1 . we arrive at 

it = —2^v — ^r]h\a\ 2 a + c.c. — |x^ 3a + c - c - 

+ ^-iuj h (7 - iuj p) a + c.c. + O (|a| 5 ) . (75b) 




1000 2000 3000 4000 5000 6000 7000 
t 

FIG. 4: An example of the kink being accelerated by the 1 : 1 
direct driving. Here h = 0.012, 7 = 1 x 10~ 3 and p = 0. 
The crosses are measured values from numerical simulations 
of the original PDE, while the lines are the predictions of the 
amplitude equations. 



B. Reduced Dynamics in Three Dimensions 

Thanks to the a-dependent driving terms in Eq. (|75bl) , 
the direct driving can sustain the translational motion 
of the kink — in contrast to the parametrically driven 
cases we have considered. FigH] shows an example of 
the kink accelerated by the 1:1 direct driving force which 
simultaneously excites the wobbling. Note that results 
from the three-dimensional system (|75l) are in excellent 
agreement with predictions of the full partial differential 
equation — not only after the dynamics have settled to 
a stationary regime but also during the transient phase. 

The detailed analysis of the three-dimensional system 
(I75p will be reported elsewhere; here, we limit ourselves 



to several basic observations. Firstly, it is straightforward 
to see that when h = 0, the trivial fixed point a = v = is 
the only attractor available in the system. Secondly, nu- 
merical simulations show that as h is increased for fixed 
7 and /J, a nontrivial fixed point bifurcates from the point 
a = v = [see FigJS] (a,b)]. For lower values of p, the bi- 
furcation is subcritical (as in the case shown in Fig[5]); for 
higher p, it is supercritical. As h approaches some crit- 
ical driving strength h c , both the \a\- and w-component 
of the nontrivial fixed point tend to infinity. Finally, in 
the region h > h Cl there are no stable fixed points. In 
this region, simulations of the system (|75p reveal a blow- 
up regime, where the functions \a(t)\,v(t) grow without 
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0.01 0.011 
h 
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FIG. 5: The hysteresis loop in the 1 : 1 directly driven cf> 4 
equation with 7 = 10 -3 and p — — 10 -4 . The driving strength 
is increased from h = 8 x 10~ 3 to 13 x 10 -3 in increments of 
0.2 x 10 -3 and then reduced back to 8 x 10~ 3 . The "crosses" 
are measured values from numerical simulations of the partial 
differential equation (|63[1 whereas the continuous and dashed 
lines show the stable and unstable fixed points of the ampli- 
tude equations (|75p . 



bounds. 

To determine the critical value h c , we assume that the 
blow-up regime is self-similar, that is, that the growth 
of v is pegged to that of a. This assumption can be 
formalised by expanding v and Arg(a) in powers of large 
lol: 



v = V 3 \a\ 3 + V 1 \a\+V^ 1 \a\- 1 - 
1-2 1 a 1 „ 1 — 4 



= oq 



-2H 



_ 4 |a|" 



. (76a) 
.(76b) 



Substituting these expansions in Eq. (|75ap and (|75b|) , and 
equating coefficients of like powers of \a\ to zero, we can 
evaluate the coefficients V n and 9 n to any order — this 
justifies the assumption. 

In particular, setting to zero the coefficients of \a\ 3 in 
Eq. (|75ap . gives 



■ C „ 2^o ICI 1 



(77) 



On the other hand, substituting Eqs.(|76p in Eq. (|75bp . we 
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FIG. 6: Results of numerical simulations of Eq. (|63[) with h 
raised beyond the critical value h c — 0.021. As in Fig[5] in 
this plot 7 = 10 -3 and p = — 10~ 4 . The main panel shows 
values of \a\ and the inset values of v as a function of h. The 
"crosses" represent measurements obtained as h is increased 
from 0.0005 to 0.05 in steps of 1.5 x 10 -3 ; the "pluses" are 
obtained as h is decreased back to 0.0005. An "asterisk" 
results when a "plus" is superimposed over the "cross" at the 
same point. 



get an equation describing the growth of |a(t)|: 
d 



dt 



\a\=x\a\ + OQa\- 1 ), 



where the growth rate 
2 



ve -ie +r) * e i6o h 



2 V 3 

Substituting for e l6 ° and V3 from (I77p . this becomes 



3?r 



:0 



2^ |CI 



h . 



(78) 



The growth of \a\ is due to the h 2 term in Eq. ([78|) which 
has a positive coefficient; the growth is damped by the 
7 term. If we reduce h keeping 7 fixed, then, at h = h c 
where 



h c , = 



8w 



ICI 2 



-,1/2 



9tt i(Cv* ~ Cv) 



7 



1/2 _ 



0.6523 7 



1/2 



(79) 



the growth rate will become equal to zero. At this point 
the blow-up regime is replaced by a stable fixed point 
- which, however, still has large values of \a\ and v. 
[Reducing h further, the fixed point will persist but the 
similarity relations (|76| will no longer be valid.] Note 
that the critical value (|79p does not depend on p. Nu- 
merical simulations of Eqs. (|75p carried out for a variety 
of 7 and p, reproduce the value of h c to high accuracy. 

As h is increased towards h c and neither \a\- nor v- 
components of the stationary point are small any longer, 
the system (|75p ceases to provide any reliable descrip- 
tion for the dynamics of the wobbler. A natural question 
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that arises here, is what dynamical regime the kink set- 
tles to for h just below h c and for h above h c . In other 
words, we want to know what happens when the the wob- 
bler is driven with the small strength h of order 7 1 / 2 [so 
that the conditions (|64|) are still in place] for which our 
finite-dimensional approximation is no longer valid. To 
answer this, we have conducted a series of numerical sim- 
ulations of Eq. (|63p with h raised from h < h c to h > h c . 
The simulations reveal that in the region inaccessible to 
our finite-dimensional approximation, the kink settles to 
wobbling with a constant amplitude which is accompa- 
nied by its translational motion with a constant velocity. 
[See FigJH]] The numerically detected values of a and v 
are of order 7 1 / 3 in this region; this accounts for the in- 
adequacy of our approximation which was based on the 
assumptions a — O (7 1 / 2 ) and v = O (7). 

The behaviour of the 1 : 1 externally driven wobbling 
kink above (and just below) the critical value is an issue 
to which we are planning to return in the near future. To 
derive the correct set of amplitude equations, we will need 
to use our asymptotic approach with modified scalings for 
a and v. 



C. 1:1 Directly Driven Wobbler and Its Radiation 

Up to O (e 2 ) , the perturbation expansion gives 

(p(x,t) = tanh [(1 - 3|a| 2 + 9/i 2 )£] + [asech£tanh£ 
+ h(l - 2 sech 2 £)] e int + c.c. + 2\a\ 2 sech 2 £ tanh£ 

- 4h(a + a*) sech£ (1 + sech 2 f) 

- /i 2 (3tanh£ + 8sech 2 £tanh£) 

+ [« 2 /i(0 + haf 3 (0 + h 2 U(0] e 2int + c.c. + (\a\ 3 ) . 

(80) 



b ± (o,o,---;7i,t 2 ,...) = -J 2 ° 



Equations (|82|) are the boundary conditions for the am- 
plitude fields B + and £>_ . The equations of motion for 

these fields are obtained at the order e 3 of the outer ex- 

(2) 

pansion. Namely, the solvability conditions for ip 3 , the 
coefficient of the second harmonic at the order e 3 , give 

D X B+ + c a dxB+ = 0, Xi > (83a) 
L>iB_ - codiB- = 0, X 1 < 0, (83b) 

where cq = ko/(2ujo). 

As before, the analysis of the linear transport equations 
(I55|) under the boundary conditions (f8"2"]l is straightfor- 
ward. The initial condition B + {X\, 0, ...) defined in 



For sufficiently large i, the variable £ is given by a; — 
vt — x , where Xq is determined by the initial conditions. 
The complex constant a and the real v are components 
of a stable fixed point (trivial or nontrivial) of the system 
(|75| . The functions fx, fz and f± are given by Eqs. (fT2|) . 
(|69|) and (l70|) . As in all previously considered driving 
regimes, the 1 : 1 direct driver excites a standing wave 
with the amplitude proportional to the driver's strength 
and frequency equal to the frequency of the driving [first 
two terms in the second line in (|8H)) ] . The standing wave 
includes also the second and zeroth harmonic, both with 
the amplitudes of order h 2 (terms in the fourth and the 
last line). 

Like the expansions in the previous sections, Eq. (|80|) 
is only valid at distances |£| = O (1). To describe the 
waveform at longer distances, we consider the outer ex- 
pansions 

4> = ±1 + e(He luJoTo + c.c.) + e 2 (f> 2 + e 3 3 + . . . 

in the regions Xq > and Xq < 0, respectively. Substi- 
tuting in the equation (|65p . the order e 2 gives 

2 = T 3iJ 2 ± ^flViwoib + c c 

+J ± B ± e l ^ ±T °- k±x ^+c.c.. (81) 

Here the top and bottom sign pertain to the Xq > 
and Xq < region, respectively. The ampli- 
tudes B± are functions of the "slow" variables: B± = 
B± (Xx , Tx , ...) and the normalisation coefficients J± 
have been introduced for later convenience. Matching 
the outer solution (|8Tj) to the inner solution ([9]) with co- 
efficients as in (|57| and and choosing J± — 
±(2 — ifco), we obtain lo± = 2wo, k± — ±ko and 



±i—k J?HA(0,0,...;T 2 ,T 3 ,...). (82) 

I 

Xi > 0, propagates, unchanged, to the right and the ini- 
tial condition B_(Xx, 0, ...) defined in Xi < 0, propa- 
gates, unchanged, to the left, both with the velocity co- 
in the expanding region —cqTx < X\ < cqTi, the ampli- 
tudes B± are constants defined by the conditions (f82|) . In 
terms of the second-harmonic radiation, this corresponds 
to two groups of radiation waves, diverging to the left 
and to the right, and leaving a sinusoidal waveform of 
constant amplitude in between. 



I 

A 2 (0,0,...;T 2 ,T 3 ,...)-^ 2 



18 



D. Qualitative Analysis 

The reason why the external force does not couple di- 
rectly to the wobbling mode in the case of the 1 : 1 
external driving (as it did in the case of the 1 : 1 para- 
metric resonance) , is the discrepancy in the parity of the 
driving profile and the wobbling mode. Instead, there 
are three indirect amplification mechanisms at work in 
this case. In each of these, the central role is played 
by the even-parity standing wave excited by the driver. 
Firstly, the square of this standing wave couples to the 
wobbling mode via the term e 3 </> 3 . Secondly, the sec- 
ond and zeroth harmonic of the standing wave as well 
as the second-harmonic radiation excited by the stand- 
ing wave couple to the wobbling mode via the term 
e 3 0o</ ) i ( / , 2- Thirdly, when the kink moves relative to the 
standing wave, the odd-parity wobbling mode acquires 
an even parity component which then couples to the 
standing wave. [This process is accounted for by the 
term e 3 FZ?o<9o0i m Eg. ([55)) .] Since the first two mech- 
anisms rely upon a quadratic superharmonic of the in- 
duced standing wave (with the amplitude of the super- 
harmonic being proportional to h 2 ), and since the veloc- 
ity of the kink (which determines the amplitude of the 
even component of the wobbling mode in the third mech- 
anism) is small, the 1 : 1 direct resonance is weak. 

VII. CONCLUDING REMARKS 

In this paper, we have used the asymptotic method 
to study the wobbling kink driven by four types of res- 
onant driving force, viz., the 1 : 1 and 2 : 1 parametric, 
and 1 : 2 and 1 : 1 external driving. We have demon- 
strated the existence of resonance (i.e. the existence 
of sustained wobbling with nondecaying amplitude de- 
spite losing energy to radiation and dissipation) in all 
four cases. This conclusion (verified in direct numeri- 
cal simulations of the corresponding partial differential 
equation) agrees with results of Quintero, Sanchez and 
Mertens who also demonstrated the existence of the res- 
onance in the 1 : 1 parametrically and 1 : 2 directly 
driven (f> 4 equation [2, 3, 4]. However, we are in disagree- 
ment with these authors on the 1 : 1 directly driven, 
damped equation. Namely, our method does capture the 
resonance in this case whereas their collective coordinate 
approach does not. (In fairness to the pioneering work 
of Quintero, Sanchez and Mertens, their direct numerical 
simulations did reveal a resonant peak at the frequency 
of the driver equal to the natural wobbling frequency of 
the kink — an experimental result which, however, did 
not reconcile with their collective coordinate predictions 
@,|.) 

In each of the four driving regimes that we have consid- 
ered in this paper, we have derived a system of equations 
for the complex amplitude of the wobbling coupled to 
the velocity of the kink. The predictions based on this 
dynamical system are in agreement with results of the 



direct numerical simulation of the full partial differential 
equation. In three out of four cases considered, the veloc- 
ity of the kink is shown to decay to zero as time advances, 
as a result of which the dimension of this dynamical sys- 
tem reduces from 3 to 2. Only in one case (the case of 
the 1 : 1 directly driven wobbler) does the velocity of the 
kink not necessarily decay to zero. In this latter case the 
wobbling of the kink is accompanied by its motion with 
nonzero velocity. 

Each of the four dynamical systems derived here give 
rise to a bifurcation diagram featuring bistability and 
hysteretical transitions in the wobbling amplitude. In the 
1 : 1 parametric and 1 : 2 direct resonances, the bista- 
bility is between two nonzero values of the wobbling am- 
plitude, whereas in the 2 : 1 parametrically and 1 : 1 di- 
rectly driven </> 4 equations, one of the two stable regimes 
involves a nonzero and the other one a zero amplitude. 
It is fitting to note here that the collective coordinate 
approach [3, H, H[ does not capture the bistability and 
hysteresis. 

In section lTVl we ranked the two parametric resonances 
according to the amplitude of the stationary wobbling re- 
sulting from the driving with a certain reference strength, 
h. Adding to this hierarchy the two direct resonances 
produces the following ranking. The 1 : 1 parametric 
resonance is the strongest of the four cases; in this case 
the amplitude of the stationary wobbling, a, is of the 
order /i 1 / 3 . The 2 : 1 parametric resonance is second 
strongest; in this case the kink responds with the wob- 
bling amplitude a ~ h 1 ^ 2 . The 1 : 2 direct resonance has 
a ~ ft 2 / 3 and the 1 : 1 direct resonance is the weakest: 
a ~ h. (The fact that the harmonic direct resonance is 
weaker than the superharmonic one, is in agreement with 
results of 0, Q where it was established in the undamped 
situation.) 

Our asymptotic approach also allows to rank the res- 
onances according to the widths of the corresponding 
Arnold tongues on the "driving strength vs driving fre- 
quency" plane. The 1 : 1 parametric resonance is the 
widest one; in this case the resonant region is bounded 
by the curve h ~ p 3 / 2 . The 2 : 1 parametric resonance is 
second widest; in this case the Arnold tongue has h ~ p. 
The 1 : 2 direct resonance has h ~ p 3 / 4 and the 1 : 1 
direct resonance is the narrowest: h ~ p 1 / 2 . We should 
also mention that the 1 : 1 parametric and the 1 : 2 di- 
rect resonance have no threshold in the strength of the 
driver whereas the 2 : 1 parametric and 1 : 1 direct reso- 
nances occur only if the driving strength exceeds a certain 
threshold value. 
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